A High Sensitivity Search for z^ e 's from the Sun and Other Sources at KamLAND 
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Data corresponding to a KamLAND detector exposure of 0.28 kton-year has been used to search 
for v e 's in the energy range 8.3 MeV < Ep < 14.8 MeV. No candidates were found for an expected 
background of 1.1±0.4 events. This result can be used to obtain a limit on i> e fluxes of any origin. 
Assuming that all v e flux has its origin in the Sun and has the characteristic 8 B solar i/ e energy 
spectrum, we obtain an upper limit of 3.7 x 10 2 cm -2 s _1 (90% C.L.) on the v e flux. We interpret 
this limit, corresponding to 2.8xl0 -4 of the Standard Solar Model 8 B v e flux, in the framework of 
spin-flavor precession and neutrino decay models. 

PACS numbers: 26.65+t,13.15+g,14.60.St,13.35.Hb 
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Of the many mechanisms that have been suggested to 
explain the solar neutrino problem jlj, neutrino oscilla- 
tions are strongly favored by the data. Assuming CPT 
invariance, the recent observation of reactor v e disappear- 
ance by the Kamioka Liquid Scintillator Anti-Neutrino 
Detector (KamLAND) |2|, combined with direct mea- 
surements of the solar neutrino flux jjj, indicates that 
the oscillation parameters lie in the Mikheyev-Smirnov- 
Wolfenstein (MSW) 4] Large Mixing Angle (LMA) re- 
gion 0. However, the limited precision of current mea- 
surements still allows for the possibility that other mech- 
anisms play a sub-dominant role. Since further study of 
the nature of neutrinos and the properties of the Sun is 
vital, we report in this Letter on a search for solar £/ e 's. 

There are several conceivable mechanisms which would 
lead to a v e component in the solar flux incident on Earth. 
Electron neutrinos with a non-zero transition magnetic 



moment can evolve into v^s or P r 's while propagating 
through intense magnetic fields in the solar core. These 
neutrinos can, in turn, evolve into 77 e 's via flavor oscilla- 
tions. There is also neutrino decay, in which a heavy neu- 
trino mass eigenstate decays into a lighter anti-neutrino 
mass eigenstate 0, ■ 

The analysis presented in this Letter concerns a search 
for DeS regardless of origin. Possible non-solar sources of 
P e 's at KamLAND include Weakly Interacting Massive 
Particle (WIMP) annihilation in the Sun and Earth || 
and relic supernova neutrinos 0, lioj . either of which 
could contribute to a continuous v e flux. The event rates 
from these an d other non-solar sources are expected 
to be small, however, and we choose to focus on models 
that predict a flux of v e 's descendant from solar neutri- 
nos. 

KamLAND was designed to study the flux of reactor 
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zVs. While the reactor v e flux spectrum has an endpoint 
of ~ 8.5 MeV, the 8 B solar neutrino flux spectrum ex- 
tends well beyond this energy to ~ 15 MeV. As a result, 
KamLAND data may be used to search for P e 's in the 
solar neutrino flux over an energy range largely free of 
reactor D e events. 

The detector consisted of a thin plastic- walled balloon, 
13 m in diameter, filled with about 1 kton of liquid scin- 
tillator (7.6xl0 31 free protons). The balloon was sur- 
rounded by an 18-meter-diameter stainless steel sphere 
instrumented with 1325 17-inch and 554 20-inch Hama- 
matsu photomultiplier tubes (PMTs), which provided 
34% photo-coverage. For the search presented here, only 
the data from the 17-inch PMTs were analyzed, lowering 
the photo-coverage to 22%. The space between the stain- 
less steel sphere and the balloon contained a mixture of 
dodecane and isoparaffin oils to act as a buffer against 
external backgrounds. The stainless steel sphere and its 
contents (hereafter referred to as the inner detector (ID)) 
was itself contained within a cylindrical water Cerenkov 
outer detector (OD) equipped with 225 20-inch PMTs. 
The OD was used to tag events due to cosmic ray in- 
duced particles. The entire detector was shielded by a 
rock overburden of about 1000m (2700m.w.e.), which 
reduced the cosmic muon flux by a factor of 10 5 relative 
to that at the surface. 

Electron anti-neutrinos were detected via the inverse 
/3-decay reaction 

v e +p — > e + + ri, (1) 

consisting of a prompt energy deposit from the positron 
and two annihilation 7's followed ^210 /is later by neu- 
tron capture on hydrogen, producing a 2.2 MeV 7. For 
the energy range of our search, this reaction, in which 
the final state neutron is free, can occur only on the free 
proton of the hydrogen nucleus. The v e energy was de- 
duced from the prompt energy E promp t using the relation- 
ship Ejj e = ^prompt +^>recoii+ 0.8 MeV, where the small 
quantity E reco u refers to the neutron kinetic energy in 
the final state and was neglected. 

Event reconstruction for high energy inverse /3-decay 
events in this analysis was similar to that described in J2J 
and was found to be accurate to within 2% from com- 
parison with the observed energy distribution of the /3 
decay of cosmogenically produced 12 B and 12 N (Fig- 
ure nj . The measured energy resolution of KamLAND 
for this data set was 7.5%/ sj E pr0 mpt (MeV) . Events with 
8.3 MeV < Ec c < 14.8 MeV, followed 0.5 /js -660 fjs 
later by a delayed event depositing between 1.8 MeV and 
2.6 MeV of energy, were selected. The distance between 
the prompt and delayed vertices was constrained to be 
less than 160 cm and both vertices were required to be 
within 550 cm of the detector center in order to suppress 
backgrounds due to natural radioactivity and muon spal- 
lation. Backgrounds were further reduced by using ID 
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FIG. 1: Prompt energy spectrum of 12 B decay. Points are 
KamLAND data and the curve is the expected f3 decay spec- 
trum convolved with the detector response. 
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solar anti-neutrino 




FIG. 2: Energy distribution of the final event candidates. 
The tail from reactor v e events is visible below 8 MeV. 

PMTs to reconstruct a muon track for all events contain- 
ing OD data. Anti-neutrino candidates associated with 
detected muons were discarded if they occurred within 

2 s after un-reconstructed muons, within 2 s after muons 
depositing at least 3 GeV, or within 2 s and less than 

3 m from a reconstructed muon track. Spallation neu- 
trons associated with tagged muons were also removed 
and did not contribute to the background. 

Figure |21 shows the delayed versus prompt energy dis- 
tribution for events after all selection cuts, except those 
on the prompt and delayed energies themselves. Taking 
into account the 12% deadtime associated with muon re- 
jection, the total sample livetime, corresponding to the 
period March 4 - December 1, 2002, was 185.5 days. 

The detection efficiency for inverse /3-decay events 
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was estimated from Monte Carlo simulation and cal- 
ibration data to be 84.2 ± 1.5%. The main contri- 
butions to the detection inefficiency were the cuts on 
the distance between the prompt and delayed vertices 
(89.8 ± 1.6%), the time between the prompt and de- 
layed vertices (95.3 ± 0.3%), neutron capture on pro- 
tons (99.5%), and the energy of the delayed event 
(98.9 ± 0.1 %). The efficiency of the vertex-separation 
cut was determined by a Monte Carlo simulation checked 
against AmBe neutron source data. The neutron capture 
time distribution with mean 210±5 /zs and the delayed 
energy cut efficiency were measured using both AmBe 
neutron and spallation neutron data. 

No events were observed in the signal region over the 
0.28 kton-year sample. A separate analysis of the Kam- 
LAND data was carried out as a cross-check using a sub- 
set of the 0.28 kton-year sample presented here. The 
results of both analyses were consistent. 

The energy spectrum of reactor P e 's extends to about 
8.5 MeV and may have constituted a small background 
in the solar v e energy region due to the detector's finite 
energy resolution. The number of background events was 
estimated to be 0.2 ± 0.2 and uncertainties in the detec- 
tor energy scale and neutrino oscillation parameters were 
the dominant sources of error. The background from 
atmospheric neutrinos was estimated, using the Barr- 
Gaisser-Stanev flux 0], to be 0.001 events. 

Cosmic ray muons interact in and near the detector 
producing spallation neutrons with an energy spectrum 
extending up to several hundred MeV. These neutrons 
constituted a source of background for the inverse (3- 
decay measurement as the prompt deposit of the neu- 
tron's kinetic energy followed ^210 ^s later by the cap- 
ture of the thermal neutron was indistinguishable from 
the inverse /3-decay event signature. As outlined below, 
we estimated the spallation neutron contribution to the 
background using a sample of neutrons selected from the 
data. 

Spallation neutron cuts were the same as for v e can- 
didates except that the fiducial volume cut was dropped 
and the muon-related cuts were replaced by the require- 
ment that at least 5 PMTs in the OD fired. The ra- 
dial distribution of the remaining candidates was fitted 
in order to obtain a smooth extrapolation of the fast 
neutrons into the fiducial volume. The resulting fitted 
function was integrated inside the volume to estimate the 
expected number of fast neutron events Nf n meeting the 
selection criteria. We used this quantity to estimate the 
two components of the fast neutron background by mul- 
tiplying Nf n by a factor of 0.11, determined from Monte 
Carlo calculations, to obtain the contribution from fast 
neutrons due to muons passing through the rock near 
the detector and by scaling Nf n by the OD detector in- 
efficiency to obtain the contribution from fast neutrons 
produced by muons passing through the OD but miss- 
ing the ID. Summed, these two components contributed 



0.3 ± 0.2 events to the background. 

We estimated the background due to accidental coin- 
cidences using data events falling within an off-time de- 
layed coincidence window of 1-10 s. Two hundred and 
seventeen such coincidences were found, corresponding to 
a background contribution of 0.02 events after normaliza- 
tion to the width of the D e delayed coincidence window. 

The residual backgrounds from cosmogenic 8 He [t\j2 = 
0.12 s) and 9 Li (t 1/2 = 0.18 s) decays were estimated by 
determining the total number of these events in the data 
sample and extrapolating into the v e signal region us- 
ing known decay times and vertex distributions. Above 
8.3 MeV, the 9 Li contribution dominated and, accord- 
ingly, analyses in that energy region dealt exclusively 
with 9 Li. The residual contribution to the background 
was calculated to be 0.6 ± 0.2 events. 

Table [I] summarizes the background estimates for this 
data set. 



Background Source 


Expected Events 


Reactor v R 


0.2 ± 0.2 


Atmospheric neutrinos 


0.001 


Fast neutrons (Nf n ) 


0.3 ± 0.2 


Accidental coincidences 


0.02 


8 He & 9 Li 


0.6 ± 0.2 


Total 


1.1 ± 0.4 



TABLE I: Estimated backgrounds for the inverse /3-decay sig- 
nal in the energy range of 8.3 MeV < E De < 14.8 MeV for 
185.5 live-days. 

The v e flux integrated over the energy range 
8.3- 14.8 MeV is obtained from: 



axexTxp p xj v 

where N S i gna i is the number of detected P e 's, 
a = 6.9 x 10~ 42 cm 2 and e = 0.84 are the cross sec- 
tion and detection efficiency, respectively, aver- 
aged over energy, T = 1.6 x 10 7 s is the livetime, and 
Pp x fv — 4.6 x 10 31 is the number of target protons in 
the fiducial volume f v (radius 550 cm). For calculat- 
ing the average cross section and detection efficiency, the 
shape of the Standard Solar Model 8 B flux without os- 
cillations Q was used. 

Systematic uncertainties in the quantities in Equa- 
tion [21 are tabulated in Table [H] The systematic un- 
certainty in the number of target protons (p p xf v ) was 
obtained by adding in quadrature the 2.1% uncertainty 
in the amount of scintillator in the balloon and the es- 
timated 3.7% uncertainty in the fiducial volume. This 
latter estimate is based on the difference between the 
measured number of spallation products in the fiducial 
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volume and the expected number assuming that the spal- 
lation products were uniformly distributed. The contri- 
bution from the energy threshold was calculated using 
the uncertainties in the energy scale (2%) and the slope 
of the neutrino flux at the threshold of 8.3 MeV. 



Quantity 


Systematic Uncertainty (%) 


Detection efficiency (e) 


1.6 


Cross section (a) 


0.2 


Number of target protons 


4.3 


Energy threshold 


4.3 


Livetime (T) 


0.07 


Total 


6.3 



TABLE II: Systematic uncertainties in quantities used to de- 
termine the flux of solar v,,. 



We derived an upper limit on 4> P( , using the 
Feldman-Cousins unified approach 15] supplemented 
with Bayesian modifications to account for the errors on 
nuisance 0] and background parameters [ItI Hsl H^ |. 
For no observed events, the upper limit of the v e flux 
was 3.7 x 10 2 cm _2 s _1 at 90% CL. Using the prescrip- 
tion described in |l5j . the sensitivity of this measure- 
ment was 7.9 x 10 2 cm~ 2 s~ 1 (90% C.L.). Normalizing to 
the solar 8 B v e flux |20| in the analysis energy window 
(8.3 MeV < E„ e < 14.8 MeV, containing 29.5% of the to- 
tal flux of 5.05xl0 6 (1.00;°; 2 g) cm" 2 s- 1 20]), this flux 
limit corresponds to an upper limit on the neutrino con- 
version probability of 2.8 xlO -4 at the 90% C.L. and rep- 
resents a factor of 30 improvement over the best previous 
measurement 21]. 

We have assumed a non-oscillatory solar z7 e flux up 
to now in order to retain as much generality as possible 
but, in the following, we have interpreted the KamLAND 
upper limit on the solar v e flux in the framework of two 
models: spin-flavor precession combined with neutrino 
oscillations and neutrino decay. 

Assuming that the solution to the solar neutrino prob- 
lem lies within the LMA region of parameter space and 
that the MSW effect is a dominant mechanism affect- 
in g th e solar neutrino flux, we followed the treatment 
of [221 |23j (taking the value of 34 degrees for the mixing 
angle) and obtained the following limit on the product 
of the neutrino transition magnetic moment \x and the 
transverse component of the magnetic field Bt in the 
Sun at a radius of 0.05i? s : 



\i B T (0.05R S 
10- 12 /iB 10 kG 



< 1.3 x lO" 3 



(3) 



The current best limit on the neutrino magnetic moment 
is from the MUNU experiment l^: fi Ve < 1.0x10 
(90% C.L.). 



Similarly, for quasi-degenerate neutrino masses, we 
were able to constrain the lifetime 0,0 for f 2 , the heavier 
neutrino, to r 2 /m 2 > 0.067 s/eV. If the neutrino mass 
spectrum is hierarchical, the limit is weaker and for m 2 
of about 0.01 eV (~ \J Ato 2 2 ), t 2 > 11 fis. Our limit, 
obtained using the appropriate decay branching ratio for 
final states containing a v e , represents an improvement 
over the current bound of r/m > 10~ 4 s/eV suggested 
in 0. 

To summarize, we have described a search for 
z/ e 's in the energy range (8.3 MeV < E 9e < 14.8 MeV) 
with KamLAND. The KamLAND detector's source- 
independent sensitivity allows for the measurement of v e 
fluxes independent of origin. No events were found in the 
185.5 live-day data set, allowing for an upper limit to be 
set on the flux from any source producing t? e 's in the ap- 
propriate energy range. We have obtained a flux limit of 
< 3.7 x 10 2 cm" 2 s _1 (90% C.L.), assuming a solar 
origin and an un-oscillated 8 B neutrino energy spectrum. 
This limit has been used to constrain models of neutrino 
spin-flavor precession and neutrino decay. 

The KamLAND experiment is supported by the COE 
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